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Abstract 

We present here the astrophysical parameters of three stellar open cluster candi- 
dates using PPMXI0 database. In this study, the main photometric, astrometry 
and statistical parameters of Ruprecht 13, Ruprecht 16 and Ruprecht 24 are es- 
timated for the first time. 
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Open star clusters (OCs) play an important role in studying the forma- 
tion and evolution of the Galactic disk and the stellar evolution as well. The 
fundamental physical parameters of OCs, e.g. distance, reddening, age, and 
metallicity are necessary for studying the Galactic disk. The strong interest of 
OCs results come from their fundamental properties. Among the 1787 currently 
OCs, more than half of them have been poorly studied or even unstudied up to 
now, Piatti et al. (2011). Thus the current paper is a part of our continuation 
series whose goal is to obtain the main astrophysical properties of previously 
unstudied OCs using modern databases. The most important thing for using 
PPMXL database lies in containing the positions, proper motions of USNO- 
Bl.cH and the Near Infrared (NIR) photometry of the Two Micron All Sky 
Survey (2MASSJ1 which let it be the powerful detection of the star clusters 
behind the hydrogen thick clouds on the Galactic plane. 

Our candidates are selected from among the unstudied OCs of Ruprecht 
catalogue (hereafter, Ru 13, Ru 16, and Ru 24). They are lying in the third 
quarter of Galactic longitude and close to the Galactic plane. The only available 
information about these clusters are the coordinates and the apparent diame- 



* Tel/Fax: +202 25560645/ +202 25548020 

1 http: //vizier, cfa. harvard. edu/viz-bin/VizieR ?-source=I/31 7 

2 http: //vizier, cfa. harvard. edu/viz-bin/VizieR ?-source=I/284 
A http: //vizier, cfa. harvard. edu/viz-bin/VizieR ?-source=II/246 



Preprint submitted to New Astronomy 



September 16, 2011 



Table 1: Equatorial, Galactic positions and estimated diameters of the clusters under investi- 
gation; sorted by right ascensions. 
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ters, which were obtained from 

webdaQ site and the last updated version of 
DIAS0 collection (version 3.0, 2010 April 30). These clusters are sorted by right 
ascensions and listed in Table 1. The quality of the data is taken into account 
and the physical properties of each cluster are estimated by applying the same 
methodology. 

This paper is organized as follows. PPMXL data extraction and preparation 
are presented in Section 2, while the data analysis and parameters estimations 
are described in Sections 3. Finally, the conclusion is devoted to Section 4. 

2. PPMXL Data Extraction and preparation 

The astrophysical parameters of the investigated cluster have been estimated 
using the PPMXL Catalogue of Roeser et al. (2010). It is combining the 
USNO-B1.0 proper motion of Monet et al. (2003) and NIR JKK S pass-band of 
2MASS Point Sources of Cutri et al. (2003). USNO-B1.0 is a very useful cata- 
logue, which gives us an opportunity to distinguish between the members and 
background/held stars. While, NIR 2MASS photometry provides J (1.25/xm), 
H (1.65/im), and K s (2.17/im) band photometry for millions of galaxies and 
nearly a half-billion stars (Carpenter, 2001). This survey has proven to be a 
powerful tool in the analysis of the structure and stellar content of open clusters 
(Bica et al. 2003, Bonatto & Bica 2003). The photometric uncertainty of the 
2MASS data is less than 0.155 at K s ~ 16.5 magnitude which is the photometric 
completeness for stars with \b\ > 25°, Skrutskie et al. (2006). 

It is noted that the candidate clusters are located near the Galactic plane 
(|6| < 10 °), therefore we expect significant foreground and background field 
contamination. These clusters have low central concentration, as appeared from 
their images on the Digitized Sky Survey 

(DSSJE 

see Fig. 1. Their apparent 
diameters are less than 10 arcmin, hence the downloaded data have been ex- 
tended to reach the field background stars, whereas the clusters dissolved there, 
i.e. the data are extracted at such sizes of about 20 arcmin. 



\http: //obswww. unige. ch/webda\ 
http://www. astro, tag. usp. br/^wilton/ 

:/ / cadewww. dao. nrc . ca/ cadebin/ getdss 



2 



Figure 1: The blue images of the candidate clusters. From left to right: Ru 13, Ru 16 & Ru 
24 respectively, as taken from Digitized Sky Surveys. North is up, east on the left. 

Comparing the data with Cutri et al.'s Point Source Catalogue (2003), we 
found that most stars' magnitudes have a "AAA" quality flag, which means the 
Signal Noise Ratio is SNR > 10, i.e. they have the highest quality measure- 
ments. Every star in our data has 3-colour photometric values J, H, K s mag; and 
proper motion values in right ascension and declination, i.e. (pm a cos S) and 
(pm 5) mas yr -1 . According to Roeser et al. (2010), the stars with proper mo- 
tion uncertainties > 4.0 mas yr _1 have been removed. Also, in this context, the 
stars with observational uncertainties > 0.20 mag are excluded, and the photo- 
metric completeness limit is applied on the photometric pass-band 2MASS data 
to avoid the over-sampling at the lower parts of the cluster's colour magnitude 
diagrams (CMDs) (cf. Bonatto et al. 2004). Pm vector point diagram (VPD) 
with distribution histogram of 5 mas yr _1 bins have been constructed as shown 
in Fig. 2. The Gaussian function fit to the central bins provides the mean pm 
in both directions. All data lie at that mean ±1 a (where a is the standard 
deviation of the mean) can be considered as astrometric probable members. In 
addition, the photometric membership criteria is adopted based on the location 
of the stars within ±0.1 mag around the ZAMS curve in the CMDs (Claria & 
Lapasset 1986). 

3. Data Analysis 

3.1. Cluster's Centre and Radial Density Profile 

The centre of any cluster can be roughly estimated by eye, but to determine 
the centre's coordinates of our candidates more precisely, we applied the star- 
count method to the whole area of each cluster. All the data that obtained 
around the adopted centre are dividing into equal sized bins in a and 8. The 
cluster centre is define as the location of maximum stellar density of the cluster's 
area. The Gaussian curve-fitting is applied to the profiles of star counts in a & 
5 respectively (cf. Tadross 2011). The estimated clusters' centres are found to 
be in agreement with that obtained from Webda with 3 s in a and 10 in 5. 
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Figure 2: Proper motion vector point diagrams VPDs for our candidate clusters after avoiding 
all data with pm errors > 4 mas yr _1 . From left to right: Ru 13, Ru 16 & Ru 24 respectively. 
Histograms of 5 mas yr _1 bins in both directions are drawn. The Gaussian function fit to the 
central bins provides the mean values of pm a cos <5 = (-2.6 ± 0.1), (-2.0 zt 0.09) & (-0.6 ± 

0. 06) mas yr" 1 and pm <5 = (4.9 ± 0.05), (5.5 ± 0.1) & (1.6 ± 0.06) mas yr _1 respectively. 

To establish the radial density profile (RDP) of a cluster, we counted the 
stars within concentric shells in equal incremental steps of r < 1 arcmin from 
the cluster centre. We repeated this process for 1 < r < 2 up to r < 20 arcmin, 

1. e. the stellar density is derived out to the preliminary radius. The stars of the 
next steps should be subtracted from the later ones, so that we obtained only 
the amount of the stars within the relevant shell's area, not a cumulative count. 
Finally, we divided the star counts in each shell to the area of that shell those 
stars belong to. The density uncertainties in each shell was calculated using 
Poisson noise statistics. 

Fig. 3 shows the RDP from the centre of the cluster to a maximum angular 
separation of 10 arcmin for the investigated clusters. To determine the structural 
parameters of the cluster more precisely, we applied the empirical King model 
(1966). The King model parameterizes the density function p(r) as: 

P(r) = h g + 1 + [° /rc)2 (1) 

where fo and r c are background, central star density and the core ra- 
dius of the cluster, respectively. From the concentration parameter c, defined 
as c = {Rum/ Rcore), Nilakshi et al. (2002) concluded that the angular size of 
the coronal region is about 6 times the core radius. Maciejewski & Niedziel- 
ski (2007) reported that Ru m may vary for individual clusters between about 
2R core and 7R core . In our case, from the c values, we can see that Ru m vary 
between 5.8R core and 6.7R core , i.e. tending to the upper limit of Maciejewski 
& Niedzielski (2007). The cluster limited radius can be defined at that radius 
which covers the entire cluster area and reaches enough stability with the back- 
ground density, i.e. the difference between the observed density profile and the 
background one is almost equal zero. It is noted that the determination of a 
cluster radius is made by the spatial coverage and uniformity of PPMXL pho- 
tometry which allows one to obtain reliable data on the projected distribution 
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Figure 3: The radial density distribution for stars in the field of the candidate clusters. 
From left to right:Ru 13, Ru 16 & Ru 24 respectively. The density shows maximum value 
at the cluster' centre and then decreases down to that point at which the decrease becomes 
asymptotically. The curved solid lines represent the fitting of King (1966) model. Errors bars 
are determined from sampling statistics (1/%/jV, where N is the number of stars used in the 
density estimation at each point). The limited radii Rn m = 4.5, 3.5 & 5.0 arcmin; and the 
core radii r c = 0.77, 0.52 & 0.83 arcmin respectively. The dark areas refer to the background 
field densities where f bg = 4.6, 6.0 & 6.8 stars per arcmin 2 respectively. 

of stars for large extensions to the clusters' halos. On the other hand, the con- 
centration parameter seems to be related to cluster age, i.e. for clusters younger 
than about 1 Gyr, it tends to increase with cluster age. Nilakshi et al. (2002) 
notes that the halos' sizes are smaller for older systems. Finally, we can say 
that open clusters appear to be larger in the near-infrared than in the optical 
data, Sharma et al. (2006). 

Knowing the cluster's total mass (Sec. 3.3), the tidal radius can be given by 
applying the equation of Jeffries et al. (2001): 

R t = 1.46 (M^ 1 / 3 (2) 

where Rt and M c are the tidal radius and total mass of the cluster respectively. 

3.2. Colour- Magnitude Diagrams 

The main photometrical parameters (age, distance, and reddening) are deter- 
mined for our clusters by fitting the 2MASS ZAMS solar metallicity isochroncs 
of Marigo et al. (2008) to both CMDs (J, J-H & K s , J-K s ) of each cluster. 
Marigo et al. (2008) isochrones can be selected and downloaded from the site 
of Padova isochrones^- Several isochrones of different ages are applied to each 
cluster; the best fit should be obtained at the same distance modulus for both 
diagrams, and the colour excesses should be obeyed Fiorucci & Munari (2003) 's 
relations for normal interstellar medium, see Fig. 4. 

The observed data has been corrected for interstellar reddening using the 
coefficients ratios 4^- = 0.276 and 4r- = 0.176, which were derived from ab- 

sorption rations of Schlegel et al. (1998), while the ratio = 0.118 was 

derived from Dutra et al. (2002). 



7 http: / /stev. oapd.inaf.it/cgi-bin/cmd 



5 




Figure 4: The net NIR CMDs of Ru 13, Ru 16 & Ru 24 for stars lying closely to the fitted 
isochrones and after removing all contaminated field stars. Age, distance modulus, E(J-H) & 
E(J-K S ) are estimated for each cluster as (1.0 Gyr, 10.8, 0.08 & 0.12 mag), (160 Myr, 12.3, 
0.22 & 0.31 mag) and (60 Myr, 11.8, 0.11 & 0.18 mag) respectively. 



Fiorucci & Munari (2003) calculated the colour excess values for 2MASS pho- 
tometric system. We ended up with the following results: e~" v = 0.309±0.130, 

= 0.485 ±0.150, where Ry = = 3.1. Also, we can de-reddened the 

distance moduli using these formulae: ■^A J -= 0.887, #^- = 0.322, then the 

° i^B — V &B-V 

distance of each cluster from the Sun (Rq) can be calculated. 

Once the cluster's distance Rq is estimated, then the distance from the 
galactic centre (Rg) and the projected distances on the galactic plane from the 
Sun (Xq & Yq) and the distance from galactic plane (Zq) can be determined, 
(for more details about the distances calculations, see Tadross 2011). 

3.3. Luminosity, Mass Functions and the total mass 

Since the central concentration of the cluster is relatively low, the determi- 
nation of the membership using the stellar RDP is difficult. Even though we 
determine the membership within a circle of radius r — 2 ~ 3 r c arcmin. By 
doing that, we obtained a more precise main-sequence in the CMDs. All of 
these stars are found very close to the main-sequence (MS) curve merged with 
some contaminated stars. These MS stars are very important in determining 
the luminosity and mass functions. 

The luminosity function (LF) gives the number of stars per luminosity inter- 
val, or in other words, the number of stars in each magnitude bin of the cluster. 
It is used to study the properties of large groups or classes of objects, such as 
the stars in clusters or the galaxies in the Local Group. In order to estimate 
the LF, we count the observed stars in terms of absolute magnitude after ap- 
plying the distance modulus as shown in Fig. 5. The magnitude bin intervals 
are selected to include a reasonable number of stars in each bin and for the best 
possible statistics of the luminosity and mass functions. From LF, we can infer 
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Figure 5: Spatial distribution of luminosity function for the clusters under investigation in 
terms of the absolute magnitude Mj. From left to right: Ru 13, Ru 16 and Ru 24 respectively. 
The colour and magnitude filters cutoffs have been applied to each cluster. The peak values 
lie at Mj = AA, 2.0 & 3.5 mag; and the total luminosity are found to be -1.9, -4.5 & 3.8 mag 
respectively. 



that the massive bright stars seem to be centrally concentrated more than the 
low masses and fainter ones (Montgomery et al. 1993). 

As known the LF and mass function (MF) are correlated to each other 
according to the Mass-luminosity relation. The accurate determination of both 
of them (LF & MF) suffers from some problems e.g. the field contamination of 
cluster members; the observed incompleteness at low-luminosity (or low-mass) 
stars; and mass segregation, which may affect even poorly populated, relatively 
young clusters (Scalo 1998). On the other hand, while, the physical properties 
and evolutions of the stars are related to their masses, so the determination of 
the initial mass function (IMF) is needed, that is an important diagnostic tool 
for studying large quantities of star clusters. IMF is an empirical relation that 
describes the mass distribution (a histogram of stellar masses) of a population 
of stars in terms of their theoretical initial mass (the mass they were formed 
with). The IMF is defined in terms of a power law as following: 

where is the number of stars of mass interval (M:M+dM) within a spec- 
ified volume of space, and a is a dimensionless exponent. The IMF for massive 
stars (> 1 Mq) has been studied and well established by Salpeter (1955), where 
a = 2.35. This form of Salpeter shows that the number of stars in each mass 
range decreases rapidly with increasing mass. It is noted that our investigated 
clusters have MF slopes ranging around Salpeter's value, as shown in Fig. 6. 

The mass of each star in the investigated clusters has been estimated from a 
polynomial equation developed from the data of the solar metallicity isochroncs 
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Figure 6: The mass function of the candidate clusters. From left to right: Ru 13, Ru 16 & Ru 
24 respectively. The slope of the initial mass function IMF is found to be T = -2.39, -2.37 & 
-2.37 respectively. 



(absolute magnitudes vs. actual masses) at a specific age of each cluster in- 
dividually. The summation of multiplying the number of stars in each bin by 
the mean mass of that bin yields the total mass of the cluster. Therefore, our 
clusters are found to have total masses of 230, 410 & 240 Mq for Ru 13, Ru 16 
& Ru 24 respectively. 

3.4- Dynamical state 

The relaxation time (T re iax) of a cluster is defined as the time in which the 
cluster needs from the very beginning to build itself and reach the stability state 
against the contraction and destruction forces, e.g. gas pressure, turbulence, 
rotation, and the magnetic field (cf. Tadross 2005). T re i ax is depending mainly 
on the number of members and the cluster diameter. To describe the dynamical 
state of a cluster, the relaxation time can be calculated in the form: 

N 

Treiax — 77~j j\7 Across j (4) 
O In 1 V 

where T cross — D/ay denotes the crossing time, N is the total number of stars in 
the investigated region of diameter D, and ay is the velocity dispersion (Binney 
& Tremaine 1998) with a typical value of 3 km s _1 (Binney & Merrificld 1987). 
Using the above formula we can estimate the dynamical relaxation time, and 
then the dynamical-evolution parameter r can be calculate for each cluster by: 

-- age (5) 



treiax 

If the cluster's age is found greater than its relaxation time, i.e. t 1.0, then 
the cluster was dynamically relaxed, and vice versa. In our case, all our clusters 
are found dynamically relaxed, where r = 5.8,6.7&6.0for Ru 13, Ru 16 & Ru 
24 respectively. 
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4. Conclusions 

Our procedure analysis has been applied for estimating the astrophysical 
parameters of three yet unstudied open clusters (Ru 13, Ru 16 & Ru 24). Hence, 
we found that these clusters have reasonable stellar density profiles, lying in 
the same absolute distance modulus, reddening range, and having IMF slopes 
around the Salpcter's (1955) value. On the other hand, the ages of these clusters 
are found to be greater than their relaxation times that infer that these clusters 
are indeed dynamically relaxed. All parameters of the investigated clusters are 
listed in Table 2. 
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